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LLLLLLLLLL III] PP 000000 LELLLLLLLL YY DDDDDDDD cece 


(om 

c 
ot 
| led 
at has 
bw bee 


LLLLLLLLLL 
LLLLLLLLLL 


ae 

ae 

et 

“wn 

nn 

“nn nn 
nn 
nn 
nn 
nn 
nn 


cr 
cr 

es 

bs 

ee nt be et fe tn ae at Bt td Oe 


-r 


15 
- Perform double floating potynented cal 16-SEP-1984 00:15:27 VAX/VMS Macro v04-00 


Edit Mistery 
Piskee ven 
LIBSPOLYD - Perform floating polynomial 


Page 


0 


K 15 
- Perform double floating polynomial cal .* ba ht 99: 


SOooooooooooooo 
OODWNAMNE WIN S OD NAN E WN 0 ODNOAU EW MOOR NOUS Win 


SWAIN AAAI NIAID PONOPINPONPONININYNY 2 2 Oe 


0 
0 
0 
0 
0 
0 
0 
0 
0 
0 


CoOoOCoooooooooao 


0:15:27 VAX/VMS Macro v04-00 Pa 1 
~SEP-1986 11:09:43 ELIBRIL.SREILIBPOLYD.MAR:1 / °° cq) 
eTITL byesrouyo - andy double floa ortng pol znontet A ert 
- IDEN File: LIBPOLYD.MAR Edit: 006 


SRAARAAAAAAA AAA AAA AEA AAAAAAAAAARAAAAAAAAAAAAARAAAAAAAAAEAEEAEEAEAAKAHTA et 


COPYRIGHT (c) 1978, 1980, 1982, 1984 BY 
DIGITAL EQUIPMENT CORPORATION, MAYNARD, MASSACHUSETTS. 
ALL RIGHTS RESERVED. 


® 
sd 
® 
® 
& 
THIS SOFTWARE IS FURNISHED UNDER A LICENSE AND MAY BE USED AND COPIED + 
Y IN ACCORDANCE WITH THE TERMS OF SUCH LICENSE AND WITH THE * 
LUSION OF THE ABOVE COPYRIGHT NOTICE. THIS SOFTWARE OR “ANY. OTHER 
IES THEREOF MAY NOT BE PROVIDED OR OTHERWISE MADE AVAILABLE TO ANY * 
R PERSON. NO TITLE TO AND OWNERSHIP OF THE SOFTWARE IS HEREBY * 
TRANSFERRED. * 
sd 

THE INFORMATION IN THIS SOFTWARE IS SUBJECT TO CHANGE WITHOUT NOTICE 
AND SHOULD NOT BE CONSTRUED AS A COMMITMENT BY DIGITAL EQUIPMENT * 
CORPORATION. . 
® 

® 

*® 

® 

® 

® 


DIGITAL ASSUMES NO RESPONSIBILITY FOR THE USE OR RELIABILITY OF ITS 
SOFTWARE ON EQUIPMENT WHICH IS NOT SUPPLIED BY DIGITAL. 
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> ABSTRACT: 
: Perform double precision floating point polynomial calculation. 
> ENVIRONMENT: User Mode, AST Reentrant 


; AUTHOR: Steven B. Lionel, CREATION DATE: 05-Oct-78 


; + 
3; FACILITY: General Utility Library 
| 
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gin 

- Put vereten number in standard format, with three digits 

in the edit snvarets JBS 16-NOV-7 

Add ‘'."' to the PSECT eireceive,. js s 3! -DEC-78 

Minor code improvements. $8 b-79 

Fix comments. SBL 31-July-1 

Use local handler to insure that exceptions other than the ones documented 
as being statuses are resignalled. SBL 25-Sept-1980 
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.SBTTL DECLARATIONS 
4 : INCLUDE FILES: 


: : SCHFDEF ; Condition handling symbols 
SSSDEF System symbols 


: EXTERNAL DECLARATIONS: 


OOooooooooooo 
DOOCOOCOOOOoOooo 
OOCCCCCOCoooo 


60 
61 -EXTRN LIBSSIG_TO_RET ; Library routine to convert 
6¢ 3 a signal to an error return 
ea SAD S Signal led conaiete 
3 = signalled con on 
0000 65 : 
0000 06 $ 
4 67 ; MACROS: 
00 68 ; 
0000 $ 
0000 0; 
0000 71 ; EQUATED SYMBOLS: 
0000 ie ; 
0000 7 
00000004 0000 74 arg = 4 3 argument 
00000008 0000 75 degree = 8 ; degree of polynomial 
0000000C 0000 76 coeff = 12 ; address of coefficient 
00000010 0000 77 result = 16 : result of polynomial 
000 78 ; table 
000 79 ; 
4 80 ; OWN STORAGE: 
000 81; 
0000 ¢ 
000 3 
it : 3 PSECT DECLARATIONS: 
00000000 56 -PSECT _LIBSCODE PIC, USR, CON, REL, LCL, SHR, - 
9000 87 EXE, RD, NOWRT. LONG 
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ave -SBTTL LIBSPOLYD - Perform floating polynomial 
: ; FUNCTIONAL DESCRIPTION: 


: 4 188°. ve provides the functionality of the VAX hardware 
nstruction POLYD to high level language users. 


The third operand points to a fe le (array) of double 
precision floating point coefficients. The coefficient of 
the highest order term of the polynomial is pointed to 

by the table address operand, i.e. the first table element. 
The table is specified with lower order coefficients stored 
at increasing addresses. 


The evaluation is carried out by Horner's method, and the 
result is stored at the location pointed to by the fourth 
operand. The result computed is: 


COOoooocooooo 


if d = degree 
and x = on 
result = CLOJ+x*(CCijJex*(Cl2]+ ... x*Cf£dJ)) 


The unsigned word degree operand specifies the highest 
numbered coefficient to participate in the evaluation. 


For further detail, refer to the VAX-11 Architecture 
Handbook for the description of POLYx. 


CALLING SEQUENCE: 


status.wlc.v = LIBSPOLYD (arg.rd.r, degree.rw.r, coeff.rd.ra, 
result.wd.r) 


INPUT PARAMETERS: 


ea 


DINE AN 9 OD NAME WIN © OD NAVE WIN 9 OD NAME WIN 9 OD NAME WO OONOUNS 


BEBE EE EMMA WIWIM RIP NINDININININID 2 SS 2 QOOTOOOOOOCOWOWONO 
PTE TE TE IETCE SOTO SOOO ICICI CIEE Ieee Ieee Ieee re rererererererererererererererererereereerereeree 


a aaa ain a ak aa ik ak i ak ka a kk en et a a  s | ss a as ss a 
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arg.rd.r - argument, ‘'x'’ in polynomial 
degree.rw.r - degree of polynomial (GEQ 0) 
coeff.rd.ra - table of coefficients, double floating 
IMPLICIT INPUTS: 
NONE 
OUTPUT PARAMETERS: 
result.wd.r - result of calculation 
IMPLICIT OUTPUTS: 
NONE 
00 FUNCTION VALUE: 
90 SS$_NORMAL - successful execution | 
0 SS$_FLTOVF - floating overflow 
Bf SS$_FLTUND = floating underflow ; 
0 SS$_ROPRAND - reserved operand, see VAX Architecture 
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D 
1of : manual for more details 
145 : SIDE EFFECTS: 
8 13 ; ALL other exceptions are signalled. 
a 
403C¢ 155 -ENTRY LIBSPOLYD, “M<IV,R2,R3,R4,R5> a Entry point, epee? int. ovf. 
4 138 3 and save R2, R3, , RS 
6D 15'AF 9E 009 138 MOVAB B“HANDLER, (FP) ; Enable local handler to process 
Boe 123 3; exceptions 
Oc BC O08 BC 04 B8C 75 $006 161 POLYD @arg(AP), - ; perform polynomial 
0000 166 a@degree(AP), - ; trap on exception to 
000D 16 acoef f (AP) 3; handler which will 
000D 164 3; unwind a return error 
000D 165 3: condition in RO to 
4444 138 3; caller of LIBSPOLYD. 
10 BC 50 7D a90p 198 mMOVa RO, @result (AP) 3: return value 
50 01 QA Baik 17 MOVZBL #1, RO 3 success status code 
04 0014 176 RET 3 return 
0015 17 
0015 174 
0015 175 HANDLER: 
0000 0015 176 «WORD 0 
0017 + #177 
0017 178 ;+ 
0017 179; If the sesieeio” is one of the documented exceptions for this routine, 
0017 180; call LIB$SIG_TO_RET to return it as a status. Otherwise, resignal. 
0017 181; #£Also, resignal Tf the depth is not zero. 
Oo? 18s 
50 08 AC DO 0017 184 MOVL CHFSL_MCHARGLST(AP), RO 3; Get mechanism vector address 
08 AO D5 O01B 185 TSTL CHFSL_MCH_DEPTH(RO) 3 Is depth zero? 
41 is ah; 186 BNEQ 90$ : If not, resignal 
51 04 aC OD 0 + f MOVL CHFSL_SIGARGLST(AP), R1 ; Get signal vector address 
50 04 Ai D0 0024 188 MOVI. CHF$L-SIG_NAME(R1), RO ; Get signalled condition 
048C 8F 50 861 0028 189 CMPW RO, #SS$_FLTOVF ; Compare conditions 
gf 13 8 D 190 BEQL 10$ : If it matches, don't resignal 
049C 8F 9 B) F 191 CMPW =sRO, #SS$_FLTUND 
é 13 0034 198 BEQL «=s«éd108 
0454 8F 0 861 036 19 CMPW RO, #SS$_ROPRAND 
i¢ (13 038 194 BEQL §~=s-«108 
04€4 8F 50. 861 3D 0=—s-«195 CMPW RO, #SS$_FLTUND_F 
08 12 te 136 BNEQ 5$ 
04 Al 049C 4 + ve 135 pov em. a5S8_FLTUND, CHFS$L_SIG_NAME(R1) ; Change fault code to trap code 
04B4 8F 8 2 04C 199 5$ CMP 8 #SS$_FLTOVF _F 
Oe 2 0051 0 BNEQ 90% 
04 Al 048C BF 3C 0053 1 MOVZWL #SS$_FLTOVF, CHFSL_S1G_NAME(R1) 
00000000'GF  6C FA 0058 06 10$ CALLG (AP); G*LIB$SIG_TO"RET ; Return signal as a status 


| IBSPOL YD 
rg form double floating pol 16 
pepe teal pi peep ae ope 
% Bes Ht MOVZUL ASSS_RESIGNAL, RO CYRIL. SRCIL1BPOLYD.mar:1 2% ¢f 
-END ; Resignal condition 
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Symbol table 66-SEP-1984 11:09:43 CLIBRTL.S REIL LIBPOLYD.MAR;1 (4) 
ARG = ° 4 

CHF SL -ACHARGLST z | 

HF SL-MCH_DEPT : 

CH SU~SIGRRGL $1 PY 4 

CHF $L_S1G_NAME 2 4 

OEF F = 

DEGREE = 

HANDLER iS R 8 

LIBSPOLYD 00 RG 

LIBS$S S16_10_RET tererere = X 0 

RESULT = 00000010 

SS$_FLTOVF = 900048 

SSS$_FLTOVE _F = 00000484 

SS$_FLTUND = 449894 

SS$_FLTUND F = 0000044 

SSS_RESIGN = 00000918 

ss$_ “ROPRAND” = 00000454 

boast eo ee ce nese meas + 
: ;_Psect synopsis ! 
PSECT name Allocation PSECT No. Attributes 
° 494544444 ' 00 ¢ 0.) NOPIC USR CON ABS LCL oe NOEXE NORD ~— NOVEC BYTE 
SABSS 0000 0°) 01 ¢ +1.) \NOPIC USR CON ABS LCL NOSHR EXE RD T NOVEC BYTE 
_LIBSCODE $0000 ( 10 02 ¢ 2.) PIC USR CON’ REL LCL SHR EXE RD NOURT NOVEC LONG 
a oo 
H ; Performance indicators ! 

Phase Page faults CPU Time Elapsed Time 

Initialization 30 00:00:00.04 00:00:02.39 
Command processing 104 00:00:00.31 00:00:01.53 
Pass i 187 00:00:02.69 00:00:11.77 
Symbol tabie sort 0 Sb SS 00:00:01.87 
Pass 2 52 00:90:00.61 00:00:01.19 

Symbol table output 4 00:00:00.02 00:00:00.03 
Psect synopsis output 3 00:00:00.01 00:00:00.01 

Cross-reference output 00:00:00.00 00:00:00.00 
Assembler run totals 381 00:00:04.09 00:00:18.79 


The working orf Limit was 1200 pages. 

21703 bytes_(43 pages) of virtual memory were used to buffer the intermediate code. 

There were 30 pages of symbol table space ah ieces d to hold 427 non-local and 3 Local symbols. 
7 source Lines were read in Pass 1, producing | object records in Pass 2. 

9 pages of virtual memory were used to define 8 macros. 
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H Macro Library statistics ! 


eee none queeeeseuvenencsesens$ 


Macro Library name Macros defined 
_$255SDUAZ8: (SYSLIBISTARLET.MLB;2 SEMEL He 
486 GETS were required to define 5 macros. 

There were no errors, warnings or information messages. 


MACRO/ENABLE=SUPPRESSION/DISABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:LIBPOLYD/OBJ=OBJ$:LIBPOLYD MSRC$:LIBPOLYD/UPDATE=(ENH$:LIBPOLYD) 
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